
function [ x,y,hx,hy ,dimX,dimY,Xp,Yp] = createMesh()
Nx1 = 10;
Nx2 = 20;
Ny  = 10;

Lx1 = 4;
Lx2 = 12;
Ly = 2;


dy = Ly / (Ny);

y(1) = dy;
hy(1) = dy;
for i=2:Ny+2
    y(i) = y(i-1) + dy;
    hy(i) = dy;
end



dx1 = Lx1 / (Nx1);

x(1) = -dx1/2;
hx(1) = dx1;

for i=2:(Nx1+2)
    x(i) = x(i-1) + dx1;
    hx(i)=dx1;
    
end

dx2 = (Lx2-max(x)) / (Nx2-.5);

for i=(Nx1+3):(Nx1+1+Nx2)
    x(i) = x(i-1) + dx2;
    hx(i) = dx2;
end


y = y - 3*dy/2;


dimX = size (x,2);
dimY = size (y,2);

for i=1:dimX
    for j=1:dimY
        Xp (i,j) = x (i);
        Yp (i,j) = y (j);
    end
end


